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SOME BASIC PROBLEMS 
SATISFACTORY general method of measuring forest 
evapotranspiration has yet to be developed. Many pro- 
cedures have been tried, but only the soil-water budget method 
and the micrometeorological methods offer any degree of suc- 
cess. This paper is a discussion of these procedures and the prob- 
lems that arise in applying them. It is designed as a reference 
for scientists and management specialists working in forest 
hydrology and meteorology. 

In forest hydrology, problems that require knowledge of the 
amount of water being lost to the atmosphere by evaporation and 
transpiration can be divided into two major groups: (1) esti- 
mating evapotranspiration from an existing fully forested sur- 
face, and (2) predicting the changes in evapotranspiration that 
would follow any given treatment. These problems can generally 
be subdivided according to the time and area scales involved. 
Water-yield problems may be of concern on an annual, seasonal, 
monthly, or even weekly basis. Flood-prediction problems may 
require knowing the changes in soil-water storage by monthly, 
weekly, or even shorter intervals. The area of the watershed 
under study can range from a few acres to thousands of square 
miles. 

At some time, apparently still in the distant future, we can 
expect that a rational procedure for estimating forest evapo- 
transpiration will be developed — based on soil, vegetation, and 
topographic parameters and soil and atmospheric variables. But 
we are still searching for basic understanding of the many proc- 
esses involved in the translocation of water from the soil to the 
atmosphere. Concurrently, we are trying to develop techniques 
that can be used to measure evapotranspiration under some limit- 
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ing conditions for application to some of the hydrologic problems 
discussed above. 

The standard device for accurate measurement of evapotran- 
spiration from agricultural crops is a weighing lysimeter. Other 
techniques are generally tested against a lysimeter. Unfortunately, 
because of the size of forest vegetation, the weighing lysimeter 
is unsuitable for forest evapotranspiration studies. In the gaged- 
watershed approach to forest hydrology, the watershed is con- 
sidered as a large unweighed lysimeter. If precipitation is meas- 
ured, and if measured streamflow represents all surface and 
subsurface water loss from the watershed, then the gaged water- 
shed gives evapotranspiration by difference for any period in 
which soil water storage at the end is the same as at the begin- 
ning. This last restriction often limits the method to annual or 
perhaps to growing-season values. Unmeasured surface and sub- 
surface losses often invalidate this method completely. 

The sapflow method, which measures the rate of water move- 
ment in the tree stem, cannot yet be used to determine transpira- 
tion rates (Swanson 1967). The quick-weighing method of meas- 
uring the weight loss of excised leaves or branches has been used 
by Rutter (1966) to obtain forest transpiration, but this method 
is laborious and gives only instantaneous values. And the method 
of enclosing a branch or whole tree in a tent and measuring the 
water vapor gained by the air is unsatisfactory because it alters the 
transpiration rate. These three methods, which are discussed by 
Swanson and Lee (1966), are subject to great sampling errors; 
and they give values of transpiration only. 

Empirical techniques have been used to estimate evapotran- 
spiration from measured evaporation loss of a small water surface 
(an evaporation pan or atmometer) or from one or two 
meteorological parameters such as air temperature, saturation 
deficit, or solar radiation. But the applicability of these or any 
other such methods proposed in the future cannot be tested un- 
less a suitable method of actually measuring forest evapotranspira- 
tion can be found. 


The application of either the soil-water budget method or the 
micrometeorological approaches to large watersheds will remain 
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impractical because of the spatial sampling involved. Evapo- 
transpiration from large watersheds must continue to be evalu- 
ated by semi-empirical techniques. However, continual improve- 
ment of these techniques must be based on studies involving the 
soil-water budget and micrometeorological methods on smaller 
areas. The soil-water budget method, in which the change of soil 
water content is measured, has been the normal method for 
experimental plot studies for a long time. 

Micrometeorological methods of determining evapotranspira- 
tion for crop, soil, and water surfaces have been developed only 
in the past 15 years. These involve measuring the energy balance 
of the surface and the turbulent transfer of heat and vapor in 
the air above the surface. The development of these methods has 
greatly increased understanding of the physical factors control- 
ling the evaporation process. 


THE SOIL-¥ 
BUDGET MET! 
The familiar water-balance equation provides one approach 


to measuring forest evapotranspiration (Burroughs and Schultz 
-1965). In its most general form for any given soil-vegetation 


volume, this equation’ can be written 


reer Ge Aw (1) 


where £ is the water output above the surface or evapotranspira- 
tion; P is the water input above the surface or precipitation; F is 
the flow in or out on the surface, or streamflow; G is the flow in 
or out below the surface, or deep seepage; and AW is the change 
in water content of the soil and vegetation volume, and of snow 
cover. The four terms on the right side of the equation must be 
measured in order to evaluate E; but under some conditions any 
one, two, or three terms may be negligible. 

We do not need to discuss the measurement of precipitation 
and surface runoff here. But note that P must be measured above 


1See Appendix for a description of the values in equations (1) to (14). 


the canopy if interception is to be included in E. Also, there may 
be surface inflow in some situations. 

The deep seepage term, G, has always been a bugaboo. Any 
loss or gain of water through the bottom or sides of the soil vol- 
ume sampled for the AW term is included. If the whole soil 
volume is relatively dry, deep seepage may be negligible; but 
movement can occur when the soil is drier than field capacity 
(Hewlett and Hibbert 1963). When a true water table is pres- 
ent, knowledge of the spatial variation of the height of the water- 
table surface and the hydraulic conductivity of the soil may be 
sufficient for estimating the inflow or outflow. Similarly, for soils 
drier than saturation, the appropriate soil-water suction gradients 
and the hydraulic conductivity will give, at least in theory, the 
inflow or outflow (Holmes 1960; Rose and Stern 1965). Prac- 
tical problems with instrumentation as well as sampling and 
theoretical difficulties, especially with hysteresis, cannot be over- 
emphasized. Several methods of indirectly estimating deep seepage 
have also been devised (Slatyer and McIlroy 1961; Willardson 
and Pope 1963). 

Much has been written about measuring soil water storage. 
Slatyer and McIlroy (1961), Taylor et al. (1961), and Cope and 
Trickett (1965) have given critical reviews of the various 
methods. For some studies, electrical resistance units may be 
satisfactory (Rutter 1964). However, present practice favors the 
neutron-scattering technique. The problem of stones in inserting 
access tubes is being overcome by use of power drills (Koshi 1966; 
Richardson 1966). Attempts to achieve satisfactory measurements 
within the surface 30 cm. (van Bavel 1961) have led to com- 
mercial surface probes and shields for depth probes (Pzver- 
point 1966), but the usefulness of surface-probe measurements in 
evapotranspiration studies has yet to be determined. 

Hewlett, et al. (1964) showed that instrument and timing 
errors in the neutron-scattering method are insignificant compared 
to the spatial sampling error for Carolina forest soils. They also 
showed that in time, repeated sampling at the same points causes 
the standard error of the change in mean water content to be 
only one-third of the standard error of the mean absolute water 
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content. This advantage is not obtained with gravimetric sampling. 

Hewlett, et al. (1964) obtained a standard error of about 0.3 
percent water content for 30 samples over several acres, each 
having 30-second timing intervals. They showed that further 
increase of the number of sampling locations or the sampling 
time would not significantly reduce the standard error. Convert- 
ing the 30 samples from percent by volume to inches of water, 
assuming a 75 cm. soil depth, yields 0.22 cm. as their standard 
error. This is of the same order of magnitude as the daily evapo- 
transpiration, so accurate daily difference in soil-water content 
appear to be impossible to detect unless the soil is extremely homo- 
geneous. The soil-water balance thus cannot give meaningful 
values of evapotranspiration for periods of less than 3 or 4 days. 

Many studies have been made of the water balance of forests, 
often to determine differences in evapotranspiration among cover 
types over weekly periods (Zahner 1955; Lull and Axley 1958). 
Sometimes inadequate attention is paid to the deep seepage 
problem, so results must be carefully evaluated. The soil-water 
balance approach provides information on rooting depth and 
soil-water content, but it gives no information about the meteoro- 
logical factors that govern evapotranspiration. 


MICROMETEORGLOGICAL 
APPROACH 


Two basic principles underlie the micrometeorological methods 
of determining evapotranspiration: the change of phase involved 
in evaporation requires an energy supply, and the removal of 
water vapor by the atmosphere is associated with a vapor-con- 
centration gradient. Two assumptions are necessary to develop 
the present theories: the evaporating surface must be considered 
two-dimensional and horizontal, with solid (or water) below 
and air above the interface; and the surface must be infinite 
in extent. These assumptions are required so that the energy and 
vapor flow is all vertical; there is no net flow in horizontal 
directions. 


Energy Balance 

The energy used in evapotranspiration, LE, is the product of 
the latent heat of vaporization, L,, and the mass loss of water, E. 
This energy must be supplied from energy sources. In fact, since 
the ideal surface is only two-dimensional, conservation of energy 
requires that as much energy reaches the surface as leaves the 
surface. The other sources and sinks of energy are the net radiant 
energy gained or lost, Rn (net radiation); heat gained or lost 
by air above the surface, H (sensible heat); heat gained or lost 
by soil and vegetation, S; and heat used in metabolism (mainly 
photosynthesis minus repiration), M. Melting or freezing of 
snow also acts as a heat sink or source, which will not be con- 
sidered here. The resulting surface energy-balance equation, when 
solved for evapotranspiration, is 


DE Ro Hees (2) 


All terms in this equation may have units of energy per unit 
surface area per unit time or, when integrated over a given time 
period, just energy per unit area. Net radiation is here positive 
for energy flow to the surface, while other terms are positive 
for flow away from the surface. For forests on a normal summer 
day, all terms are then inherently positive. But in other conditions, 
any one or more of the terms can become negative: L,E during 
dew formation, Rn at night, H at night or during advection (see 
below), S at night or in the autumn, and M at night due to 
respiration. 


Table 1.—Diurnal energy balance of a young fir fores# 


Flux Time Daily 

1-2 a.m. 7820.12] pm. — Se6ipan, Total 

Calorie cm-* min-} Calorie cm-* 

Rn == ().06 O55 Ce 0.38 586 
L,E =! 44 .74 22 386 
H — O02 .06 38 ao. 197 
S —s.02 .00 .03 .0O1 3 
M come Od .O5 ‘Oz —=—=~.03 0 


1Baumgartner (1956). 


Table 1 shows the values of these energy-balance components 
for a fir forest on a clear summer day (Baumgartner 1956). 
Other types of forest, and, indeed, other agricultural surfaces, 
show energy balances similar to this. The S and M terms are 
generally small for closed forests; and they may be measured or 
estimated with sufficient accuracy, or else ignored. 

The net radiation is actually a combination of two types of 
radiation — solar and thermal—each of which is transferred 
both toward and away from the surface (downward and upward). 
Since most of the energy required for evaporation is supplied by 
solar radiation, R,, evaporation is more closely related to R, than 
to any other single variable. The fraction of R, that is reflected 
by the surface is called the albedo, a. Thermal radiation, also 
called terrestrial, infrared, or longwave radiation, is emitted both 
downward by the atmosphere, R,,, and—in usually larger amounts 
—upward from the surface, R,,. The upward component is pri- 
marily emission from the surface, which is governed by the sur- 
face temperature, T,, but also includes a reflected part of the 
downward thermal radiation. Combining these radiation fluxes 
suitably gives the net radiation 


Pe Cl — ee eR eo Tt (3) 


where « is the absorptivity or emissivity of the surface for thermal 
radiation and o is the Stefan-Boltzmann constant. 

Net radiation can be measured directly with a suitable radiom- 
eter; but the other two important terms in the energy balance, 
L,E and H, are difficult to measure directly for forests. These 
two are closely related to the turbulent transfer processes in the 
atmosphere. Thus it is not possible to find evaporation by the 
energy balance alone; we must consider the aerodynamic prop- 
erties of the atmosphere. 


Aerodynamic Equations 

Four important quantities are transferred vertically by move- 
ment of the air above any surface: heat, water vapor, momen- 
tum, and carbon dioxide. In each case, the flux in question is 
proportional to the vertical gradient of, respectively, temperature, 
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water vapor concentration, wind speed, and CO, concentration. 
The CO, flux is important in plant growth studies but does not 
concern us here. 

Sensible heat is heat transferred to or away from the surface 
because the surface is colder or warmer than the air above. Mo- 
mentum transfer from the air to the surface is due to friction 
exerted on the air moving over the surface and is evidenced by 
lower wind speed closer to the surface. Water vapor generally 
moves away from the wet surface into the drier air as evapora- 
tion; but sometimes it moves in the reverse direction, resulting 
in condensation as dew or frost. 

The equations governing these transfer processes are similar:? 


H = — G, pK, (dT / dz) (4) 
E = — pK, (dq / dz) (5) 
rinpitien (dug ele) (6) 


The sensible heat flux, H, and the water vapor flux, E, appeared 
above in the energy-balance equation (2). The momentum flux 
is t. The density of the air, p, and its specific heat, C,, are con- 
sidered constants. The transfer coefficients or diffusivities, K,, 
K,,, and K,, for heat, water vapor, and momentum, respectively, 
are not constant but vary with windspeed, surface roughness, 
height, and sensible heat flux, H. The terms in parentheses repre- 
sent the vertical gradients or change with height z, of air tempera- 


w? 


ture, T; specific humidity, q; and wind speed, u. 

The three transfer coefficients are closely related, since all are 
dependent on the turbulent properties of the atmosphere. Many 
experiments have shown that 


br ane ku, 7 Sf om (7) 


where k is the von Karman constant, ¢,, is a stability correction, 


2The gradients here and later in this report can be written as total derivatives 
(dT/dz) rather than as partial derivatives (8T/5z) because horizontal homogen- 
eity and steady state conditions are assumed. When this assumption does not hold, 
advection occurs (see page 00). And to be strictly correct, the adiabatic lapse 
rate, 0.01°C m”, should be subtracted from dT/dz. 
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and u, is a form of the momentum flux called the friction velocity, 
which is defined as 


Wa (rfp)? (8) 


When the sensible heat flux, H, is zero (neutral or adiabatic 
condition), the stability parameter, ¢, is unity. This case is the 
regime of fully forced convection in which turbulence is pro- 
duced only by friction due to wind movement over the rough 
surface. When H is positive (unstable or lapse condition), free 
convection due to rising warm air and sinking cold air increases 
the turbulence, which effectively decreases ¢ to some value less 
than unity. When H is negative (stable or inversion condition), 
the turbulence is damped and ¢ increases to a value greater than 
unity. 

By an analogy called the similarity hypothesis, the other two 
diffusivities are similarly defined as 


fe keuSz, / pn 
Rae Kutz / ow 


Many scientists assume ¢ = ¢, — ¢, — ¢,, and thusK—K, 
— K, — K, as a satisfactory first approximation. This is certainly 
true in neutral conditions when all ¢ values are unity. But the 
question of equality of the transfer coefficients in diabatic condi- 
tions is an unsolved problem in micrometeorology. Several in- 
vestigations indicate that K, may be twice as large as K,, in ex- 
tremely unstable conditions and only half as large in extremely 
stable conditions (Drimmel 1968; Record and Cramer 1966). 
Researchers had thought that K, was approximately equal to K,, 
(Deacon and Swinbank 1958), but recent investigations by Dyer 
(1967) indicate that K, — K,. Until differences are quantita- 
tively established, equality will have to be assumed for practical 
work. 

A second unsolved micrometeorological problem is the form 
of ¢ under varying stability conditions. Usually ¢ is given as a 
function of one of four other stability parameters: the Richard- 
son number in either its flux or gradient forms (R, and Ri), and 
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the Monin-Obukhov dimensionless height in either its flux or 
gradient forms (z / L and z / L’). These are defined here for 
reference: 


Ree oe SN ee 8 ee 
Ty Cant (day dz) 1K — K,, T (du / dz)? 
z = ek) KK, Zien came 1. 1/taZ) 


—_— 


LemiCe pune iy Ki, Liga eid /idz) 


where g is the acceleration of gravity and T is the air tempera- 
ture. The terms are further interrelated since 


Ron — 2700 


All the parameters R,, Ri, z / L, and z / L’ increase in absolute 
value nearly linearly with height. The Richardson number at 
any height z can be found by: 


Ri= | oT ee!) In a z 


T (ug — u,)? 21 


with; I¢theimean ofS candela 

When the transfer coefficient is to be determined from wind 
and temperature measurements, the form of ¢ must be known. 
Although theory requires functions of ¢ involving R, or L, these 
cannot be used in practice unless H is independently evaluated 
through use of such techniques as the eddy correlation method 
(see below). Therefore, practical forms of ¢ are functions of 
Ri and L’ and thus involve implicit assumption of K,, = K,. 
Monin and Obukhov (1954) introduced a linear relation between 
¢ and z / L’. However, this relation contained an empirical co- 
efficient that has been found to vary with stability (Pandolfo 
1963; O’Brien 1965), so a number of more complex forms have 
been proposed (Panofsky 1963; Swinbank 1964; Webb 1965). 
Brutsaert (1965) has given an excellent summary of these relations 
and their use in evaporation measurements. Bernstein (1966) and 
others who have used existing data to compare these forms find 
that the experimental data is not good enough to show that one 
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form is better than another. Thus, in practice, nearly any of the 
forms may be chosen. In cases of extreme stability and instability, 
no proposed forms will fit the data adequately. 
The friction velocity, u,, can be found from the integrated 
form of the combined equations 6, 7, and 8: 
k (u, — uy) 


 BRONGEL aCe ear - 


where y = / aie, ee 
ors 


and u, and u, are wind speeds measured at two heights above 
the canopy, z, and z,. ¥ is the integrated form of the stability cor- 
rection, which depends on z, and z,. The zero-plane displacement, 
_D, is a correction to the measured heights of the sensors above the 
ground, z, and z,. D is about equal to the general height of the 
vegetation and allows for the fact that the top of the canopy, and 
not the ground level, influences the wind above the trees. 

If we extrapolate this wind profile toward the surface, a, 
approaches zero as z, — D approaches a constant, z,. The par- 
ameter z,, is called the roughness length and is characteristic of 
the physical roughness or turbulence-inducing properties of the 
surface. If z, and D are known for forests, u, can be obtained by 
measuring u at only one height. D and z, must be found from 
measuring wind speed at several heights, and very accurate meas- 
urements are required (Robinson 1961; Covey 1962; and Pan- 
ofsky 1963). For forests z,, is of the order of 1 m., while the 
zero plane displacement is somewhat less than the height of the 
trees ( Baumgartner 1956; Rauner 1960). 


Mass-Transfer or Profile Method 


Substituting (9) into (7) and then into (5) and integrating 
gives the mass-transfer equation for E as 


B= — ie (uz — u,) (q2 — qi) (10) 


[ = (2=B)—- +] 


sa 


This equation can be used directly for finding evapotranspira- 
tion. It requires wind and vapor concentration measurements at 
two heights above the surface and temperature measurements 
for the stability correction. The measurement accuracies required 
are high, since the differences are small. Baumgartner (1956) 
found that this method worked well over a forest at night, but 
it was unsatisfactory during the day because of insufficient meas- 
urement accuracy. Deacon and Swinbank (1958) have proposed 
a modification that eliminates the temperature measurements for 
the stability correction and assumes a constant but undetermined 
roughness length. 

Bulk aerodynamic methods use surface values, u, — 0, (2, — 
D), == z,,and.q,=—4,, in the form 


(qe a qo) (11) 


f (uz) (q2 — qo) 


This method involves only wind speed at one height but requires 
constancy of surface roughness. The exacting measurement of z, 
must be made at the outset, or else the wind function f (u,) can 
be obtained empirically. The stability correction, if not ignored, 
requires measurement of temperature at two heights. The surface 
humidity, q,, cannot be measured directly except over water, where 
it is equal to the saturation specific humidity at temperature T,. 
We should note here that 


0.622 


E 


where p is atmospheric pressure and e is vapor pressure in the 
same units as p. 


Eddy Correlation Methods 

Vertical transport of heat and water vapor in the atmosphere 
is turbulent in nature. The air flow is made up of a series of 
eddies, which are carried past a fixed point by the horizontal 
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wind. The vertical component of the wind speed, w’, at a point 
will vary continuously as the eddies pass, alternately being directed 
upward and downward. When evaporation is occurring, the air 
moving upward in the eddies will be wetter than the air moving 
downward. If we define q’ as the instantaneous departure of q 
from its mean value at the point, then the product w’q’ is the 
deviation of the instantaneous water-vapor flow from its mean 
value. The evaporation is found by averaging over time as shown 
by the overbar 


E= pwd’ 


An instrument called the Evapotron has been developed to find 
E by this method (Dyer and Maher 1965). It shows great promise 
as a research tool but has not yet been tried over forests. 

The sensible heat flux is similarly dependent on turbulent flow, 
but here temperature takes the place of humidity. 


H= C, pw’ T’ 


The Fluxatron (Dyer et al. 1957) measures H in this way but also 
has not yet been applied to forests. 


Difference Method 

The energy balance equation (2) alone can be used to 
determine E only if the sensible heat flux is found independently; 
then H is subtracted from measured Rn to give E. I have called 
this the difference method, as it apparently has no other com- 
monly used name. 

The required sensible heat flux can be found by the eddy cor- 
relation method, or it can be found by using equations analogous 
to (10) and (11) for E. Substituting (9) into (7) and then into 
(4), and integrating, gives the profile equation for H as 


= eee (u, —u,) (T, — T,) (12) 
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The difficult measurement of (q2 — q,) required for (10) its 
substituted for by the easier measurements of (T, — T,) and 
Rn. Accuracy is high over surfaces where evaporation 1s ma and 
H is small. 

The bulk aerodynamic equation for sensible heat flux is sim- 
tldietov Chee 


(Tz a @ (13) 
2 


Combining (13) with (2) is the bulk-difference method for 
evapotranspiration. A modification of this method has been used 
over a hardwood forest in Russia (Rauner 1960; Dzerdzeevski 
1963). However, no independent checks of the results were made. 
The surface temperature, T,, may be measured directly by in- 
frared radiation a (Fuchs and Tanner 1966). ony 
ments applied to equation (11) also apply here. 

Tanner (19606) has proposed a modification of (12) similar 
to Deacon and Swinbank’s (1958) modification of (10). 


Bowen Ratio Method 


The Bowen ratio is defined as H / L,E. This ratio can be ob- 
tained from (10) and (12) and substituted into (2) to give, 
after rearrangement, 


‘Ete Rn — S—M 
Gt) 
es eee 


Use of the Bowen ratio (the term in brackets) eliminates any 
need for wind measurements, known heights, stability corrections, 
roughness length, or zero-plane displacement; but the difficult 
vapor concentration difference is required (Suomi and Tanner 
1958; Fritschen 1966). The ratio K, / K, should also appear in 
the bracketed term if it is not assumed equal to unity. The Bowen 
ratio method has also been called the energy balance method, but 
this name is misleading. Baumgartner, (1956) in his classic work 
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on a fir forest in Austria, used this method, but no independent 
check on the evaporation data was obtained. 


Combination Methods 

Potential evapotranspiration is defined as the evaporation from 
a surface for which q, is the saturation specific humidity at 
T,. With this assumption, suitable combination of equations (2), 
(11), and (13) leads to 


Bee (A/G) (Rn —S —M) + ph (Gro — q2) (14) 
P+(AL,/G) 


eats 


ey 


where qr. is the saturation specific humidity at T, and A is the 
slope of the saturation specific humidity-temperature curve (van 
Bavel 1966; Tanner and Fuchs 1968). With an empirical wind 
function for h, this is the well-known Penman (1948) method. 

The “generalized combination method’ for actual evaporation 
developed by Tanner and Fuchs (1968) is simply the bulk dif- 
ference method described above as the combination of (13) and 
€2): 

A commonly used method of estimating E is to determine E, 
by (14) and then to use an empirical relationship of E / E, to 
available soil water (Denmead and Shaw 1962; Crawford and 
Linsley 1963). As soon as this relationship is put on a theoretical 
basis (Visser 1963; Cowan 1965; Woo et al. 1966), the combina- 
tion method for potential evapotranspiration and soil-water meas- 
urement may give an accurate measurement of evapotranspiration. 

The combination method and soil-water approach are also the 
avenues to predicting effects of forest treatments on evapotrans- 
piration, since they include all relevant variables. Micrometeoro- 
logical studies of forests have already led to analysis of differences 
in evapotranspiration from different cover types (Baumgartner 
1967; Rauner 1965). All the relevant surface variables — albedo, 
surface temperature, roughness, stomatal resistance, root and soil 
water properties — as well as all necessary meteorological vari- 
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ables — will some day be included in a complete quantitative 
description of the whole soil-plant-atmosphere system. 


PRACTICAL DIFFICULTIES 


instrumentation 

The measurements needed for any of the micrometeorological 
methods unfortunately require sensitive instruments that are not 
as yet standardized. Net radiation can be measured with a variety 
of different instruments, some of which are commercially avail- 
able (Gates 1965; Federer 1967). Differences in wind speed, 
temperature, and vapor concentration between two heights require 
extremely accurate sensors because the differences are small. The 
required order of accuracy needed for the differences are 1 cm. per 
sec. for wind, and 0.02° C. for wet-bulb and dry-bulb tempera- 
ture. Greater height differences allow lower absolute accuracy in 
the measurements. Anemometers meeting these requirements can 
be purchased. Temperature and water vapor sensors of sufficient 
accuracy are usually built by the investigator, but no commonly 
used design is available. An outstanding review of such instru- 
mentation is given by Tanner (1963). 

If methods utilizing surface values of temperature and vapor 
concentration are used, the required accuracies are of a magnitude 
less than those given above. However, these methods require 
either an empirical wind function or determination of roughness 
length and zero-plane displacement from accurate wind profiles, 
and the assumption that these remain constant. 

All the methods discussed here require time integration of the 
meteorological variables. Integration can be done either elec- 
tronically at the site before recording, or by digital computer on 
a series of recorded, instantaneous values. Integration times may 
vary for different studies from 15 minutes up to half a day. De- 
velopment of satisfactory recording systems is a significant elec- 
tronics problem, especially for forest measurements where remote 
operating stations that may be unattended for periods of a week 
or more may be required. Calibration drift in either the sensors 
or the electronics is a serious problem for the systems that must 
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be employed. Frequent interchange of sensors between heights 
is one way to reduce this problem (Tanner 196042). Compromises 
between the drift problem and unattended operation may be 
needed. 

In short, a great deal of work is needed on instrumentation 
problems before micrometeorological methods can be easily ap- 
plied to forests. 


Advection 


The current theory of the micrometeorological evaporation 
methods assumes a horizontally uniform, two-dimensional surface. 
If this assumption is not met, there will be a net horizontal trans- 
fer of energy into or out of the volume being measured. This hori- 
zontal transfer is called advection, and it is not accounted for by 
the theory. The three-dimensional form of forests and their hori- 
zontal variation imply that advection may be a serious problem in 
applying micrometeorological methods. 

Absorption of the horizontal component of solar radiation is 
exemplified by isolated trees or stands that have a sunny side and 
a shadow. Hilly topography also produces inhomogeneities in hor- 
izontal solar and infrared transfers. Geiger (1965, pp. 369-381) 
and Lee (1964) have analyzed these situations. Lateral intercep- 
tion of water droplets as fog drip is an extreme example of hori- 
zontal advection of water vapor, but it can also occur in invisible 
forms. 

Advection of sensible heat and water vapor occurs in three 
forms (fig. 1). Net movement into or out of a stand, or of in- 
dividual trees below the tree tops, is called the clothesline effect 
(Tanner 1957). There is no theory for this effect, and its magni- 
tude has never been measured. Local clothesline advection pro- 
duced by individual tree crowns and canopy openings causes local 
variation in energy balance and in mass transfer. Spatial integra- 
tion then becomes necessary and can be achieved either by using 
several sensors at different locations or by raising the sensors 
higher above the canopy. 

When sensors are elevated to obtain better spatial sampling 
and to obtain maximum possible vertical differences, a second 
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Figure 1.—The three kinds of advection, or net lateral 
energy transfer. 


form of advection becomes important. Energy can be transferred 
laterally in the space between the top of the canopy and the top 
of the highest sensing instrument. With increasing height, the 
wind, temperature, and vapor concentration are affected more and 
more by surfaces farther upwind and less by the surface imme- 
diately below the sensor. The measured vertical gradients are then 
not representative of the underlying surface. The fetch require- 
ment of a site is the minimum upwind distance at which a sig- 
nificantly different surface may exist without affecting conditions 
at the sensor height. This type of lateral transfer may thus be 
called fetch advection. 

Several studies have been made of fetch requirements in the 
case of a dry surface upwind of a wet surface (de Vries 1959; 
Philip 1959; Dyer 1963; Rider et al. 1963). Fetch advection can 
also be produced by a change in surface roughness that affects 
the downwind turbulent transfer coefficients (Elliott 1958; Brooks 
1961; Panofsky and Townsend 1964). A rough rule of thumb 
requires a fetch 50 to 100 times the highest sensor height above 


18 


the zero plane. However, no measurements are available to test 
this rule over forests. The forest meteorologist is seriously 
squeezed between the minimum height allowed by clothesline ad- 
vection and the maximum height allowed by fetch advection. 

Energy can also be transferred to an area by horizontal transfer 
above the sensor height, then by vertical transfer to the surface. 
This third type of advection is called the oasis effect because the 
condition is most pronounced over desert oases. Oasis advection 
is satisfactorily measured by micrometeorological methods since 
the transfer is vertical at sensor height. 

Clothesline advection is net horizontal transfer below the tops 
of the trees. Fetch advection is net horizontal transfer between 
the tree tops and the sensing instrument height. Oasis advection 
is net horizontal transfer above sensor height, or above treetop 
height if there are no sensors. The theoretical equations discussed 
in this paper are correct only for regions of height in which none 
of these three advection forms is appreciable; that is, only for 
regions in which all net transfer is vertical. 

A form of vertical heat transfer within the volume being con- 
sidered must also be mentioned. Studies by Denmead (1964) 
and Knoerr and Gay (1965) indicate that radiant energy absorbed 
at the upper surface of the forest canopy is in part converted to 
sensible heat, which is transferred downward into the canopy and 
used in increasing evaporation. The effective source level, and 
thus the zero-plane displacements for sensible heat, latent heat, 
and momentum, may be different. The effect of this phenomenon 
on the various micrometeorological evaporation methods has not 
yet been studied. 


CONCLUSIONS 


Forest hydrology research and forest watershed management 
have been hindered by lack of a satisfactory method for measur- 
ing forest evapotranspiration. Soil-water sampling used in the 
water balance is limited to periods greater than 3 days and is 
plagued by the problem of unmeasured deep seepage. Micro- 
meterological methods, based on the energy balance and vapor 
transport, offer an approach to measuring evapotranspiration for 
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periods as short as 15 minutes. However, only two successful 
studies applying these methods to forests have been reported 
(Baumgartner 1956; Rauner 1960). 

Several different micrometeorological methods can be used to 
obtain nearly independent results. Future experiments should 
evaluate the usefulness and accuracy of these methods. The magni- 
tude of advection error due to surface inhomogeneities, and con- 
sequent spatial sampling requirements, need study. Theoretical 
problems involving equality of the eddy diffusivities and the dia- 
batic correction may not be serious for'the accuracy required here. 
But problems with instrumentation require further research. 

Combining micrometeorological theory with results from soil 
physics and plant physiology should lead to methods of predict- 
ing the effects of changes in the soil-plant-atmosphere system on 
evaporation. 
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APPENDIX 


The units given below represent a consistent set, but they are by no 
means the only units that are commonly used in the soil-water budget 
method and the micrometeorological methods. Energy-flux values are more 
commonly given with time units in minutes rather than seconds, as shown 


a albedo dimensionless 
C, specific heat of air (= 0.24) cal g? °K? 
D zero plane displacement cm 
E evapotranspiration (rate) 2 cl sec 
e evapotranspiration (quantity ) cm 
[note £ = eee E dt] 
Pw 

E, potential evapotranspiration ecu so. 
e vapor pressure mb 
F streamflow cm 
G deep seepage cm 
g acceleration of gravity (= 980) cm sec? 
H _ sensible heat flux cal cm-? sec 
h wind function cm sect 
K diffusivity cm? sec-1 
k von Karman constant (= 0.42) dimensionless 
L = Monin-Obukhov flux length cm 
L’ Monin-Obukhov gradient length cm 
L, heat of vaporization (= 580) cal o* 
M heat for metabolism cal cm=? sec 
P precipitation cm 
p atmospheric pressure mb 
q specific humidity 

(g water / g air) dimensionless 
q’ instantaneous deviation from 

mean specific humidity dimensionless 
Gr2 saturation specific humidity 

at temperature T, dimensionless 
R, flux Richardson number dimensionless 
Ri _—_— Richardson number dimensionless 
R,q downward thermal radiation cal cm-? sec" 
Ri, upward thermal radiation calm sec 
Rn net radiation Caron sec 
R, downward solar radiation cal cm-? sec"? 
S change in soil or vegetation 

heat storage cal cm-? sec"? 
ag temperature Og 
T’ instantaneous deviation from Ox 

mean air temperature 
u—s mean wind speed cm sect 
u, _ friction velocity cm sect 


od 


w’ __ instantaneous vertical 


component of wind cm sec-} 
z height cm 
z, ‘roughness length cm 
A _ slope of saturation specific 
humidity-temperature curve Ox 
AW change in soil water content cm. 
é thermal emissivity dimensionless 
p air density : (1.0012) 9 an 
pw water density (21.00) 9 cm 
o Stefan Boltzmann constant (=.5.7:\X 10°) 2 sec 
= momentum flux Bictn sec 
¢@  diabatic correction _ dimensionless 
y integrated diabatic correction dimensionless 


Subscript , refers to the surface; , and , refer to heights above the 
surface. 

Subscripts ,,, », and , on ¢ and K refer to momentum, heat, and 
water vapor, respectively. 
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THE FOREST SERVICE of the U. S. Depart- 
ment of Agriculture is dedicated to the principle of 
multiple use management of the Nation's forest re- 
sources for sustained yields of wood, water, forage, 
wildlife, and recreation. Through forestry research, 
cooperation with the States and private forest 
owners, and management of the National Forests 
and National Grasslands, it strives —as directed 
by Congress—to provide crea gienter 
service to a growing Nation. 


